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Abstract. In the 1990s a comparison of sparse EGRET measurements with single-dish flux density monitoring from 
the Metsahovi and UMRAO programs established a temporal connection between the onset of flaring at radio band 
and the occurrence of gamma-ray activity. Correlations between the emergence of new VLBI components from the 
core, flares in linearly polarized radio flux, and gamma-ray activity in bright EGRET-detected blazars supported 
a picture in which the gamma-ray and the radio band emission arises in the same shocked region of the jet, with 
the high energy emission produced via inverse Compton scattering by the synchrotron-emitting electrons in the 
jet. Quantitative tests of this scenario, however, were hampered by insufficient temporal sampling of the data and 
the simple nature of the models adopted. The extensive data from Fermi coupled with the wealth of well-sampled 
radio band data from old as well as new programs such as the F-GAMMA project now permit statistical studies for 
large numbers of sources, including weak HBLs, and detailed analyses of individual highly-active class members. 
I summarize progress in understanding the origin of the gamma-ray emission using these new measurements. I 
focus on three areas: attempts to isolate the physical site of the high energy emission using time delay information; 
investigation of the emission process using the characteristics of the variability; and quantitative tests of the shock 
model picture using high-time-sampled multifrequency linear polarization data, VLBP imaging, and new models 
of propagating oblique relativistic shocks incorporating detailed radiative transfer calculations. 



1. Introduction 

In spite of the limitation provided by the relative low sen- 
sitivity and poor time sampling (detection statistics only) 
from EGRET, these measurements set the framework for 
our current understanding of the MeV-to-GeV 7-ray emis- 
sion from AGN. The enormous increase in sensitivity and 
the use of large-area sky scan mode rather than point- 
ings made possible by the > 2 sr field of view accessi- 
ble to the Large Area Telescope (LAT), the principle in- 
strument for studying blazars onboard Fermi has result- 
ing in 7-ray light curves with sufficient time sampling to 
track the variations on intraday time scales during ex- 
treme strong flares, and on daily-to- weekly time scales for 
tens of bright AGN in the energy band 20 MeV to > 300 
GeV. To match these data, a new generation of radio band 
monitoring programs has emerged, providing well-sampled 
light curves for hundreds of sources. These enable investi- 
gators to carry out both statistical studies of large samples 
of AGN to low radio-band flux limit, and detailed studies 
of individual objects using multifrequency data provided 
in support of the Fermi mission. Using well-established 
procedures, including cross-correlations of infra-band data 
permitting direct comparison of the activity across bands 
and investigations of the character of the variability itself 
through analyses which probe the emission process, de- 
tailed studies of the variability properties in both bands 
can be carried out. These comparison rely on single dish 
monitoring data which provide temporal resolution com- 
plementing the spatial resolution provided by monthly-to- 
less-frequent VLBI imaging measurements. I summarize 



the current status of these monitoring programs and re- 
view what they can tell us about the site of and mechanism 
for the generation of 7-ray emission from AGN detected 
by Fermi. 



2. The EGRET era: setting the stage 

The launch of the Compton Gamma Ray Observatory 
(CGRO) in April 1991 provided us with a new view of fa- 
miliar blazars. One of the early discoveries from EGRET 
(the instrument onboard operating above 30 MeV) was 
the identification of a population of GeV-7-ray-emitting 
blazars, positionally associated with well-known radio- 
bright AGN. Comparison of activity monitored in the mil- 
limeter band (37 & 22 GHz) in a sample of « 70 bright 
sources combined with the EGRET detections in the first 
EGRET catalogue (Fichtel et al. 1994), the results of an all 
sky survey carried out from 1991 April to November 1992, 
established a statistical connection between enhanced ra- 
dio band activity and detection in the energy band E>100 
MeV (e.g. Valtaoja & Terasranta ll995[) . Note that only 
a handful of photons were detected by EGRET above 
10 GeV (Thompson I2006p . and that much information 
was extracted based on only a small number of pho- 
tons in these early studies. Even with the paucity of 
EGRET data and the relatively small number of sources 
detected, it became quickly apparent that these ^/-ray- 
detected sources were predominantly radio-bright, com- 
pact, core-dominated blazars with high brightness temper- 
atures (Valtaoja & Terasranta 1996j). Further, a temporal 
association between the occurrence of an EGRET dctec- 
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Radio flare phase 
(0=start, 1=maximum, 2=end) 

Fig. 1. Time of occurrence of 7-ray detection as a function Fig. 2. 7-ray flux versus Metsahovi peak flare flux. Figure 
of flare phase based on Metsahovi data. Figure courtesy courtesy A. Lahteenmaki. 
A. Lahteenmaki 



tion and the rise portion of the radio band flare (start of 
flare to maximum) was established for individual events 
using both the millimeter (Valtaoja & Terasranta 1996) 
and the centimeter band (Aller, Aller & Hughes |1996[) 
total flux density monitoring data. An example plot il- 
lustrating this result is shown in Figure 1 which includes 
EGRET phase 1+2+3 7-ray data and Metsahovi monitor- 
ing measurements obtained during this time period (see 
also Lahteenmaki & Valtaoja 1999). Because of the severe 
undersampling of the EGRET data, however, it was not 
possible to carry out cross-correlation analyses between 
the two bands. 

Attempts to correlate 7-ray photon flux and radio flare 
peak flux or radio luminosity yielded mixed results that 
ranged from no clear correlation (Aller, Aller, & Hughes 
1996) to claimed correlations. This discrepancy stemmed 
from a variety of causes including the use of non simulta- 
neous data, dependence of the luminosity on redshift, and 
time delays between the radio and 7-ray activity produced 
by opacity in the radio jet (Miicke et al. 11996)) . An example 
flux-flux correlation is shown in Figure 2 based on detec- 
tions from EGRET and monitoring observations from the 
Metsahovi program (see Lahteenmaki et al. H997[) . 

Studies of a handful of individual, exceptionally-bright 
sources for which the binned EGRET measurements 
yielded multiple detections (time-averaged photon flux 
values as a function of time), in combination with radio 
band light curves and VLBI imaging data, confirmed the 
physical association between flaring in the two bands iden- 
tified statistically. Analysis of the data for sources such 
as NRAO 530, confirmed the general trends but found 
cases of 7-ray detections during declining or constant 
millimeter-to-entimeter band flux (Bower et al. 1997) 
which deviated from the general trends. 



The link to flaring, and subsequently the identification 
of a specific flare phase established a connection to the de- 
velopment of individual radio band outbursts. Since the 
radio band variations are generally attributed to shocks 
in the radio jet, this connection led investigators to the 
conclusion that shocks might also play a role in the gen- 
eration of the 7-ray flaring (e.g. Valtaoja & Terasranta 
IT996)) . 

Linear polarization monitoring data track tempo- 
ral changes in magnetic field structure. In well-resolved 
events, such data can be used to search for the signature of 
a shock: a swing in electric vector position angle (hereafter 
EVPA) and an increase in the fractional linear polariza- 
tion (hereafter LP). These changes are associated with the 
compression of the initially tangled magnetic field during 
the passage of the shock. The 22 & 43 GHz VLBA imaging 
data for 42 7-ray-bright blazars obtained from the Jorstad 
& Marscher program, in combination with UMRAO lin- 
ear polarization monitoring data, were used to identify a 
correlation between high 7-ray state and an increase in 
polarized radio band flux supporting a picture in which 
the radio emission and the 7-rays originate in the same 
shocked region of the jet, with the 7-ray emission pro- 
duced by inverse Compton scattering in the thin forward 
region of the shock where the electrons are accelerated 
(Jorstad 200T]). This scenario placed the 7-ray emission 
site parsecs far from the central engine and downstream 
of the radio band core. 

The majority of the EGRET-detected sources were 
highly luminous HPQs. With only a few exceptions, these 
were the very same radio-bright sources which had been 
included for decades in the Metsahovi and UMRAO flux 
monitoring programs. We now know that this was a se- 
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lection effect resulting from the limited sensitivity of 
EGRET. 



3. Single dish radio band monitoring during the 

Fermi era 

In spite of the progress made in understanding the origin 
of the 7-ray emission using the EGRET data, answers to 
many questions remain in areas where single dish moni- 
toring data can provide insights. Specific questions which 
can be addressed using these data are: 

1. Where within the jet is the 7-ray flare produced? This 
question can be examined using cross identification of 
flares in the high energy and radio bands. A combination 
of single dish monitoring and VLBI imaging can be used 
to localize the physical site of the high energy emission 
using the radio core as a fiducial reference point. 

2. What emission mechanism is responsible for the gener- 
ation of the 7-ray emission? This can be explored using 
SEDs based on contemporaneous broadband measures, 
including single-dish total flux density data at multiple 
epochs. These trace the evolution in the vF v versus v 
plane. 

3. What is the mechanism for acceleration of the parti- 
cles? Tests for the presence of shocks during 7-ray flaring 
using a combination of data and modeling can be carried 
out. 

4. What special conditions are present in the jet during 
broadband flaring? Information on jet properties during 
individual flares and on changes in jet conditions can be 
derived using single dish measurements. 

Fermi scans the sky every 3 hours providing data 
which can be analyzed to study the variability in the 
MeV-GeV band on near-daily to weekly time scales. In 
support of this program several new programs as well as 
continuing old ones monitor the total flux density of AGN. 
The salient features of several of these programs provid- 
ing data in the centimeter-to-millimeter band for sources 
north of approximately -30° are summarized in Table 1. 
The first three entries form the F- GAMMA alliance which 
works in tight coordination. The OVRO program provides 
time-sampling at a single frequency, 15 GHz, temporally 
matched to that attained using Fermi. Observations in 
this band are well-suited for tracking variability since the 
amplitude of the variability is relatively large compared to 
that at lower frequencies, and the measurements are rel- 
atively free from tropospheric effects which increase with 
higher frequency. The UMRAO and Metsahovi programs 
provide long term histories for smaller numbers of sources 
and with less-frequent cadence. The AGN population now 
detected by Fermi includes a substantial population of 
radio-weak sources belonging to the HBL class. These 
radio-weak sources can only be monitored using single 
dish instruments with large collecting areas such as the 
OVRO 40-m telescope and the Effelsberg 100-m instru- 
ment. Spectral data is provided by Effelsberg/IRAM and 
the RATAN-600 programs. Together the data base of ob- 
servations includes both temporal and spectral coverage 



Table 1. Current Single Dish Monitoring Programs 



Prog. 


freq. (GHz) 


sampling 


size/advantage 


OVRO 


15 


2-3/week 


many 


>1000 
sources; low S 


Effsbg. 


2.64-43 


monthly 




^60 
spectra 


IRAM 


86-270 


monthly 




^60 
inner jet 


UMRAO 


4.8, 8, 14.5 


1-2/week 


35 

mf; 


in core grp. 
includes LP 


Metsahovi 


37 


monthly 




«100 

inner jet 


RATAN-600 


1-22 


2-4/year 




600 
spectra 



for a large number of sources. Multifrequency data, includ- 
ing linear polarization measurements, are themselves im- 
portant for identifying both self- absorption and Faraday 
effects. 

4. Is the broad band activity causally related? 
Radio band-7-ray correlations 

Evidence in support of a statistically-significant correla- 
tion between the fluxes in the two bands would support the 
hypothesis that the emissions are related. Results based 
on observations from the OVRO program for 49 sources 
in the Fermi LAT three month catalogue are shown in 
Figure 3 from the work of Richards et al. (|2009|) . The data 
shown are averages over the three month period August 
4 - October 30, 2008. A correlation is apparent visually. 
A formal analysis gives a correlation coefficient of r=0.56; 
Monte Carlo simulations indicate a chance probability of 
this correlation of only 5xl0 -4 . 

Figure 4 shows the flux-flux correlation obtained using 
results from MOJAVE for sources in the high-confidence 
3-month LAT list (Abdo et al. 121305)) with radio flux den- 
sity >0.2 Jy and S >-30 degrees. The 15 GHz fluxes shown 
include both MOJAVE data (total flux density over the 
entire VLBA image) and supplementary UMRAO and 
RATAN-600 fluxes (shown by open circles). This work 
found a correlation between the 7-ray flux in the first three 
months of Fermi operation and these quasi simultaneous 
15 GHz fluxes, and identified that the radio jet was in an 
enhanced state near the time of strong 7-ray emission for 
sources in this group (Kovalev et al. [2009). 

5. Characterization of the Variability 

Variability time scales, variability amplitude, quasi- 
periodicities, and indicators of the noise process producing 
the variability can be extracted from the detailed study 
of light curves. Common analysis tools include structure 
functions, PDSs, and wavelets. Useful measures of the de- 
gree of variability include the fluctuation index and the 
normalized excess variance. Similar characteristics in the 
variability of the two bands would support a related ori- 
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Fig. 3. Fermi LAT flux density in the 100 MeV energy 
band versus OVRO 15 GHz flux density based on data 
obtained in the first 3 months of operation (Richards et 
al. 12009)1 . 
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Fig. 4. Average 7-ray photon flux in the 100 McV-1 
GeV energy band versus 15 GHz flux for sources in the 
MOJAVE sample. 



gin. However, a sufficiently long data train is required to 
capture the characteristics of the variability. 

5.1. Results 

Structure functions, autocorrelations and power spectra 
are related measures of the distribution of power versus 
time. Structure functions can provide both a characteristic 
time scale and a "measure" of the noise process; the for- 



mer comes from the time lag where the curve plateaus and 
the latter from the slope of the curve. First order struc- 
ture functions were computed for the best observed (i.e. 
most variable) sources in the UMRAO program (Hughes, 
Aller, & Aller I199"2"|) : the turnover yielded a characteristic 
variability time scale of 2 years while the slope, gener- 
ally near 1, was consistent with shot noise. Because of the 
length of the data train in the radio band, several events 
in each source were captured in general; in a few sources, 
e.g. 3C 84 and 3C 454.3, the variability was dominated by 
a single event, and in these cases the slope of the structure 
function was 1.55 -1.6. 

In the case of the 7-ray data, nearly two years of data 
have accumulated. A study containing light curves of 106 
sources showing weekly averages of the data (E>300MeV) 
obtained over the first 11 months of operation has recently 
been submitted for publication which includes analysis us- 
ing structure functions and PDSs. (Abdo et al. I2010c|) . 
Characteristic structure function power indices are in the 
range 1.0 to 1.6. Higher values are found for 0235+164 
and 3C 454.3; in these sources the light curves are dom- 
inated by a single large event. The light curves of these 
two objects in the radio band show the same characteristic 
behavior during this time period. 



5.2. Time-Delay Studies 

The site of the 7-ray emission can in principle be iden- 
tified relative to the emission site for the radio emission 
from time delays by using well-sampled light curves to 
look for correlated activity. If similar features can be un- 
ambiguously identified, the associated time delays, provide 
information on the location of the 7-ray emission site. In 
the radio band this can be referenced to the radio core 
at 43 GHz identified from VLBI imaging. Monitoring ob- 
servations from the OVRO program have traced the total 
flux density variability at 15 GHz in hundreds of sources 
with a time sampling well-matched to that provided by the 
LAT. Cross-correlation of the OVRO 15 GHz and early 
Fermi data for a few representative cases are presented in 
Max-Moerbeck [20091 These few examples illustrate that a 
range of behavior patterns are present in the time delays. 
While activity is clearly associated across the spectrum, 
unambiguously identifying individual events across the 
spectral bands is extremely difficult except during rare, 
large-amplitude, resolved flares. In other cases, frequency- 
dependent factors including self-absorption and opacity, 
differ from source to source and affect the time delay re- 
sults. An underlying problem in such comparisons is that 
persistent trends may not exist from source to source, or 
even in the same source from epoch to epoch; this lack of 
persistence has recently been quantitatively examined in 
the optical to 7-ray regime and attributed to variations 
in a number of parameters which are important in the 
emission process (Bottcher & Dermer feOlOp . 
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Fig. 5. Weekly averages of the total flux density for OJ 
287 from the UMRAO program showing triple-frequency 
data obtained over 4 decades. 



5.3. Periodicity 

Quasi-periodicities have been identified in several sources 
observed in the UMRAO program using a variety of anal- 
ysis tools. The strongest evidence for periodicity in the 
centimeter band for any source is for OJ 287. Visual in- 
spection of the light curve for this source shown in Figure 5 
reveals a long term trend with superimposed flares. Quasi 
periodic behavior has been identified in this source us- 
ing wavelets (Hughes, Aller, & Aller 119981 ). The total nux 
density amplitude is currently approaching the maximum 
value attained over the UMRAO multi-decade time pe- 
riod, and the source is also bright at 7-ray band (see 
Figure 8). No evidence for periodicity has been found (yet) 
using the Fermi light curves over the first 11 months of 
operation, and it will be very interesting to see what anal- 
yses of longer data trains reveal. 



6. Emission Process: studies of SEDs using single 
dish data 

With the high time sampling now available, the evolu- 
tion of the SED during individual events can be followed; 
such data can provide insight into the emission mechanism 
using modeling. The low frequency hump is generally at- 
tributed to synchrotron sclf-Compton emission produced 
in the jet, while the second (high frequency) hump is often 
attributed to a complex mix of components (e.g. Abdo et 
al. I2010dp . An example SED during flaring is shown for 
3C 279 (Abdo et al. l2010a|) in Figure 6. In the radio band 
the SEDs at the two epochs overlay, while significant dif- 
ferences in the two time intervals are apparent in all other 
bands. 
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Fig. 6. Evolution of the SED during 7-ray flaring. The 
filled symbols denote data obtained during the first five 
days of a sharp 7-ray flare. The unfilled triangles show 
data obtained about 70 days later during a time when the 
7-ray flux was relatively quiescent. In the radio band, the 
total flux density showed a small systematic decline. The 
figure was kindly provided by M. Hayashida. 
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Fig. 7. EVPA swing detected by MOJAVE in 1502+106 
at epochs bracketing the 2008 August flare. The most sig- 
nificant changes occurred during the time interval between 
the first and the second images shown. Figure provided by 
Y. Y. Kovalev. 



7. Mechanism for particle acceleration: Shocks in 
the radio jet during 7-ray flaring 

Independent evidence in support of a scenario in which 
shocks propagating outward in the jet flows are the ori- 
gin of the observed blazar activity in the radio through 
optical spectral regions has come from both model fits of 
broadband spectral data (e.g. Marscher & Gear 11983]) and 
of fits to centimeter band multifrequency flux and linear 
polarization data during outburst using transverse shock 
models (e.g. Hughes, Aller, & Aller 1198T)]) . One possibility 
for the origin of these disturbances is Kelvin-Hclmholtz- 
induced instabilities. Both detailed hydrodynamical sim- 
ulations of jet flows and stability analyses have confirmed 
that such instabilities develop naturally and commonly 
in the kinetically-dominated portion of the jet; current 
driven instabilities have also been considered and may 
dominate in the magnetically dominated regime near the 
central engine (e.g. Hardee 1200pp . Shocks have been cited 
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Fig. 8. Shock signature during flaring in OJ 287 in late 
2009. From bottom to top weekly-averaged 7-ray pho- 
ton flux and daily averages of radio band total flux den- 
sity and LP (fractional linear polarization and EVPA), 
For comparison VLBA fluxes at 43 GHz (denoted by 
stars) and a 15 GHz from MOJAVE (source-integrated 
values filled squares; core only unfilled). 7-ray fluxes 
(units photons/sec/cm 2 xlO -7 ) were kindly provided by 
S. Jorstad (private communication). 



as a means of particle acceleration resulting in flaring in 
the 7-ray band, but searches for specific evidence of shocks 
during 7-ray flaring have not been widely conducted. In 
principle, shocks in the parsec scale regime of the jet can 
be identified by searches for the expected swing in the 
orientation of the EVPA, apparent in consecutive VLBP 
images, but in general the VLBI sampling is too sparse to 
'detect' these transitions. An example of this kind of flip in 
EVPA is shown in Figure 7 which displays the MOJAVE 
15 GHz images bracketing flaring detected by Fermi 'near' 
to a major flare which occurred in 2008 August (Abdo 
et al. [2010b). In general the VLBA imaging data is too 
sparsely sampled to permit following a shock event in 
detail. Single dish monitoring of the linear polarization 
can provide the temporal resolution, but the shock signa- 
ture will be smoothed out if there are contributions to the 
emission from many competing, evolving components. In 
sources dominated by a single or only a few components, 
the signature of a shock appears in the single dish monitor- 
ing data as as ordered swing in EVPA (through 90 degrees 
in the special case of a transverse shock) combined with 
an increase in the degree of linear polarization. In the gen- 
eral case, shocks are expected to be oriented obliquely to 
the flow direction. The emission properties of such shocks 
have not been explored in detail in past studies. 

As part of the Michigan monitoring program, we are 
currently obtaining single dish monitoring observations 



Fig. 9. Simulated light curves from radiative transfer cal- 
culations showing from top to bottom total intensity, frac- 
tional polarization, and EVPA for three frequencies. The 
Lorentz factor of the flow=2.5, the Lorentz factor of the 
shocks 6.7, and the viewing angle=10 degrees. Other 
model parameters are given in the text. 



with a cadence of 1-2 times weekly of both the total 
flux density and the linear polarization for a small sam- 
ple of sources known to be flaring in the 7-ray band 
and sufficiently bright (S>0.4 Jy) in the radio band to 
look for a shock signature during times of 7-ray flaring. 
The program sources (3C 66A, 0235+164, 3C 84, 0420- 
014, 0454-234, 0528+134, 0716+714, 0727-115, 0805-077, 
OJ 287, 3C 273, 3C 279, 1308+326, 1502+106, 1510- 
089, 1633+382, OT 081, 2022-077, BL Lac, CTA 102, & 
3C 454.3) are all members of MOJAVE, and 15 of the 
21 are in the Boston U. 43 GHz blazar program led by 
A. Marscher and S. Jorstad. These VLBI data provide 
crucial complementary information on flow direction rela- 
tive to the magnetic field orientation and a pattern speed 
characterizing the jet flow, as well as a means for identify- 
ing sources with relatively simple structure. Example light 
curves obtained for one of our program sources, OJ 287, 
are shown in Figure 8. The LP exhibits the expected 
oblique shock signature - a swing in EVPA through about 
40° and an increase in fractional linear polarization from 
near to a maximum value approaching 8% - tempo- 
rally associated with the 7-ray flaring in late 2009. The 
LP swings in EVPA that we have found in this and other 
sources take place over a time period of order one to two 
months and require frequent measurements to properly 
track the variability. In addition to the UMRAO light 
curves, we show in this figure MOJAVE 15 GHz and BU 
43 GHz measurements for comparison. These illustrate the 
excellent agreement between the single dish and VLBA 
results and graphically show the sampling rate relative to 
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the rapidity of the event itself. Note that in spite of the 
self-absorption apparent in the UMRAO light curve in the 
4.8-14.5 GHz band, that the spectrum is relatively trans- 
parent in the 15-43 GHz domain as indicated by the nearly 
equal flux amplitudes in total flux and linear polarization. 

To explore the effect of the various input parameters on 
the predicted light curves and eventually to use the mod- 
els to analyze specific events, we have developed new ra- 
diative transfer models which incorporate oblique shocks. 
These models assume an extreme relativistic equation of 
state and are determined principally by two free parame- 
ters: the shock compression and the direction of the shock 
(forward or reverse). The direction is important for time 
delay considerations. From modeling of such events, in- 
formation can be extracted about the conditions in the 
jet during 7-ray flaring; changes in these conditions from 
event to event can be identified and explored as more data 
accumulates. 

We show in Figure 9 representative, simulated light 
curves from this work. In the example shown, a shock ori- 
ented at an obliquity of 45 degrees to the flow direction 
has been introduced into the relativistic flow at t=0. A 
forward-moving shock and a compression of 0.7 have been 
assumed. The model predicts the general features appar- 
ent in the data (see Figure 8): a total flux outburst, an 
increase in the linear polarization to near 10%, a swing in 
EVPA through about 40 degrees, and the general spectral 
behavior. Figure 10 contains a series of images showing the 
evolution of the simulated structure. Images such as these 
will be generated for comparison with VLBP imaging data 
to serve as model constraints in the fitting procedure. 

8. Future Work 

Comparison of the variability detected in the Fermi and 
radio spectral bands continues to support a generic sce- 
nario in which the same disturbance is responsible for both 
the radio band flares and the 7-ray activity, most likely 
a shock, located near the radio core in the spine of the 
radio jet (e.g. A. Pushkarev, This Workshop). However, 
some nagging exceptions to this simple picture remain. 
For example it does not account for the presence of rapid 
(hourly) variability in the 7-ray emission detected during 
strong flares by Fermi (Tavecchio et al. I2010|) . and known 
from earlier EGRET results. This rapid variability implies 
a small source size and can most easily be explained by an 
emission site near the central engine. Such results call into 
question the completeness of a picture which attempts to 
explain all emission as originating at one site and by one 
mechanism. Exceptions to the rule need to be reconciled 
with the generally-adopted picture. 

The specific conditions in the jet leading to 7-ray flar- 
ing have not been isolated. These can be examined using 
a combination of the unique temporal resolution provided 
by single dish monitoring, spatial information from VLBI 
studies, and detailed modeling of the emission patterns. 
Model development, such as the shock modeling described 
here, is only in its very early stages. Refinement will in- 
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Fig. 10. A sequence of simulated images from radiative 
transfer calculations showing polarization contours and 
EVPA orientation as a function of time. 



elude the development of more sophisticated models and 
fits of these models to specific flares. 

Data from large telescopes such as the OVRO 40-m 
and the Effelsberg 100-m dishes probe the variability prop- 
erties of pure BL Lacs and the large population of radio- 
weak HBLs, with sub- weekly cadences which can not be 
provided by VLBI imaging alone. Detailed analyses of 
these sources may reveal other properties and behavior 
patterns which need to be included in the general picture. 
Recent work from Metsahovi suggests that even the RBLs 
may have different emission properties characterized by 
7-ray photon fluxes which are weak relative to those from 
other source types and uncorrelated with radio-band flar- 
ing (Leon-Tavares, This Workshop). 

Long-term light curves in both the optical and ra- 
dio spectral bands, have identified that there are changes 
in the character of the variability on time scales of sev- 
eral years to decades. Searches for associated changes in 
the character of the radio band and 7-ray activity as the 
Fermi data accumulates would be further support of an 
intimate connection between the emission in these two 
regimes. Such studies will rely on historical data from 
long-term programs which provide consistently acquired 
data at a common observing frequency in combination 
with the data acquired during the operation of Fermi. 

The first indications of an association between the 
emission in the radio and 7-ray bands came from observa- 
tions of a few extreme sources and limited photon statis- 
tics by EGRET. This general picture can be refined and 
filled in using the wealth of data now available from both 
the single dish data in the radio band and Fermi data. 
This will not an easy task because of the complexity of the 
emission patterns and the dependence on many parame- 
ters, but exciting, new results will surely result enhancing 
our understanding of Blazars during the Fermi era. 
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